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ABSTRACT 

We demonstrate the possibility of detecting tidal stripping of dark matter subhalos 
within galaxy groups using weak gravitational lensing. We have run ray-tracing simu- 
lations on galaxy catalogues from the Millennium Simulation to generate mock shape 
catalogues. The ray-tracing catalogues assume a halo model for galaxies and groups, 
using various models with different distributions of mass between galaxy and group 
halos to simulate different stages of group evolution. Using these mock catalogues, 
we forecast the lensing signals that will be detected around galaxy groups and satel- 
lite galaxies, as well as test two different methods for isolating the satellites' lensing 
signals. A key challenge is to determine the accuracy to which group centres can be 
identified. We show that with current and ongoing surveys, it will possible to detect 
stripping in groups of mass 10 12 M -1O 15 M Q . 
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1 INTRODUCTION 

In the standard picture of hierarchical structure formation, 
larger dark matter halos are built up through the accretion, 
stripping, and mergers of smaller halos. At the extremes of 
the halo mass spectrum, namely isolated field galaxies and 
galaxy clusters, we have a relatively good picture of how the 
mass within these structures is organized. For isolated galax- 
ies, most of the mass is contained within a halo of dark mat- 
ter, as confirmed by galaxy-galaxy weak gravitational lens- 
ing measurements jBrainerd et al. | ll996l: [Hudson et al.lll998j: 
Smith et al1l200ll; iGuzik fc SeliaMl200a iMandelbaum et al] 



20061 : iGavazzi et ail 120071 ). Simulations have shown that 



the shape of this hal o can be well-appr o ximated by an 
NFW density profile (|Navarro et all Il997l ; iLu et all 120061 1 
and galaxy-galaxy weak gravitatio nal lensing measurements 
have confirmed these r esults l|Kleinheinrich et al.l 120031 : 
iMandelbaum et alj |2008h In galaxy clusters, most of the 
mass also seems to lie within an NFW dark matter halo, with 
the constituent galaxies co ntributing only small perturba - 
tions to the density profile (|Mandelbaum et alj|2006l J2008). 
Gravitational lensing measurements have confirmed that the 
halos around individual galaxies within clusters are signifi- 
cantly smaller than the halos around comparably-luminous 
field galaxies, and this effect is more extreme with galax- 
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ies closer to the cent res of clusters (|Limousin et al.l 120071 : 
iNataraian et alj|2009h . 

However, between the extremes of field galaxies and 
rich clusters, the picture is less clear. Since multiple galax- 
ies must merge together to eventually form clusters, at some 
point the mass in the galaxies' individual halos must migrate 
into a shared halo. This process most likely occurs through 
tidal stripping: when two galaxies pass near each other, the 
halo of the less massive galaxy will tend to be "stripped" 
from it and thus join the more massive galaxy's halo. This 
effect has been demonstrated in various N-body dark mat- 
ter simulations (iHayashi et al ] |2004l :l Kazantzidis et al 112004 



ISpringel et alj [2008h Tidal stripping is also expected to re- 
move hot gas from less massive galaxies, which will have 
the effect of quenching th eir star formation in a process 
known as "strangulation" jBalogh fc Morris! l200Ch . Galax- 
ies in dense environments are known to be significantly 
redder on average than field galaxies |Butcher fc Oemlen 
1 19841 : iMoore et al.ll 19961 : iBalogh et alJll999l . l2004 h and tidal 
stripping may contribute to the quenching of star forma- 
tion, so there is a strong motivation t o understand the me- 
chan i cs and timing of tidal strip ping |van den Bosch et al.l 
l200St iKawata fc Mul chacv 20q1). It remains unclear, how- 
ever, whether this process is rapid or gradual. 

In this paper, we focus on galaxy groups, an intermedi- 
ate mass scale between field galaxies and clusters (typically 
structures in the mass range (10 12 Mq <J Mhaio <J 10 14 M@ 
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are considered groups, and more massive structures are con- 
sidered clusters). Weak gravitational lensing provides the 
only practical tool to measure the shapes and masses of dark 
matter halos around indi vidual galaxies within groups. Lens- 
ing a n alyses of groups dHoekstra et al. 2001 ; Parker et aU 



2005 : iMandelbaum et all 120061: 



Johnston et al 



20071; 



Hamana et al.l 20091 ; Leauthaud et al. 20ld ; Ford et all 
20121 ) have shown that the group lensing signal can be 



measured and is consistent with an NFW density profile. 
However, as we show in Section 14.21 the signal could result 
from either a single, smooth halo, or a model which includes 
both a central halo and an individual halo for each satellite. 
As such, it is necessary to measure at the lensing signals 
around satellites themselves to get a full picture of the 
mass distribution. Only limited wo rk has been done in the 
group regime to date. For example, ISuvu fc Halkolal (|2010h 
studied a strong-lensing system and determined that tidal 
stripping did seem to occur around the satellite studied, 
which lies in a group of mass on the order of 10 12 M Q . 
While this result is promising, a broader base of data will 
be needed to develop firm conclusions. 

Upcoming data fro m surveys such as the GAMA-II sur- 
vey |Driver et alj|201lh will provide an unprecedented op- 
portunity for us to use weak gravitational lensing to ana- 
lyze a large sample of galaxy groups. In anticipation of the 
availability of this data, in this paper we investigate various 
methods for detecting the extent of tidal stripping of satellite 
halos. We hope to determine the best analysis method and 
to generate a framework through which we can interpret fu- 
ture results, and assess the effects of var ious sources of bias 
and error. Previous theoretical studies llMoller et al] 120021 ; 
iNataraian et al] |2002| ; lYang et al] l200rJ) have looked into 
this issue. In particular. IPastor Mira et al.l (|201 ll ) demon- 
strated a method - called "Galaxy-Galaxy Lensing with Cal- 
ibration" (GGLC) by them and the "Mirror Method" by us 
(see Section 13.3. 1[) - that can be used to isolate the lensing 
signals of satellite galaxies. Our work expands on these past 
results by testing and comparing different methods for iso- 
lating satellites' lensing signals from the contamination of 
group halos, assessing the effects of errors in estimating the 
position and mass of the group's halo on these results. To do 
so, we simulate the effects of weak gravitational lensing on a 
sample of mock galaxy catalogues extracted from the Mille- 
nium Simulation, using models that correspond to different 
amounts of tidal stripping. We measure the lensing signal 
around both satellites and group centres, and apply two dif- 
ferent methods to isolate the lensing signal of the satellites' 
halos by removing the estimated signal contribution of the 
groups' halos. By looking at the resultant lensing signals for 
satellites of similar mass, but varying the stripping model, 
the mass of the groups that contain these satellites, and 
the magnitudes of sources of error, we can determine which 
method is superior, and in which regimes it might be possi- 
ble to distinguish between stripping models. 

It is possible to study galaxies in a group environment 
using either photometric or spectroscopic redshift estimates. 
Since photometric redshift estimates are significantly worse 
than spectroscopic estima tes, with typical errors of ~ 0.05 
in z for 5-band surveyrs (|Bem'tezl I2OO0I ; iHildebrandt et al] 
2012), it is thus much more difficult to generate a group 
catalogue of galaxies. Various methods exist to determine 
the positions of groups using only photometric redshifts 



jKoester et al.ll2007l; iLi fe Yee 2008; Milkeraitis et alj|20ld : 
iGillis fc Hudsonll201ll ~ but it is difficult to determine which 
galaxies belong to thes e groups and which are field galaxies 
(Gillis fc Hudsonll201ll ). Spectroscopic redshifts provide the 
advantage that group membership can ty pically be deter- 
mine d to ~ 80% purity and completeness (|Robotham et al] 
2011). Although it is possible to investigate the presence or 
absence of tidal stripping using only photometric redshift es- 
timates, these methods are beyond the scope of this paper, 
which deals only with the case where spectroscopic redshifts 
are available. Methods based only on photometric redshifts 
will be discussed in a future paper. 

In Section [2] of this paper, we discuss the galaxy cat- 
alogues we use and our models for halo mass distribution. 
In Section [3] we explain how we simulate lensing signals 
and how we remove the influence of groups' signals from 
the lensing signals detected around satellites. In Section 21 
we present the results of our simulations and discuss their 
implications. 

For consistency with the Millenium Simulation, we use 
the following cosmological parameters: Ho = 73 km s _1 
Mpc" 1 , fi m = 0.25, 0.x = 0.75, and D. h = 0.045. All stated 
magnitudes are in the AB system. Since there is no clear di- 
vision between galaxy groups and galaxy clusters, we use the 
terminology "galaxy groups" throughout this paper, even 
though some of the structures we refer to as such would be 
more commonly deemed clusters. When masses are quoted 
in this paper, M is used to refer to the halo mass of a group, 
and m is used to refer to the halo mass of a satellite, unless 
otherwise specified. When radial measurements are used in 
this paper, R refers to a projected, 2D distance, and r refers 
to a 3D distance. 



2 SIMULATIONS AND MODELS 

In this section we discuss the simulations used to generate 
our galaxy catalogues and our models for the mass distri- 
butions of galaxies within these catalogues. In Section 12.11 
we discuss how we extracted our galaxy catalogues from 
the Millennium Simulation and the properties of these cat- 
alogues. In Section ^. 2l we discuss the mass distributions we 
assume for dark matter halos surrounding the galaxies in our 
catalogues; Section [2.2.11 discusses the mathematical form of 
the model we use, and Section 12.2.21 discusses the different 
models we use for the presence or absence of tidal stripping, 
along with other datasets we use for comparison purposes. 



2.1 Simulations 

In order to assess the differences varying stripping models 
have between their lensing signals, it is easiest to work with 
a simulated galaxy catalogue, where exact redshifts, masses, 
and mass distributions can be known. In this paper, we used 
a semi-analy tic galaxy cata l ogue based on the Millennium 
Simulation i|Springel et al] l2005h by [5e Lucia fc Blaizotj 
l|2007h . The Millennium Simulation is a collisionless sim- 
ulation of N = 2160 3 dark matter particles, each with a 
mass of 8.6 x 10 s /i _1 M , in a box with 500 h 1 Mpc sides 
and periodic boundary conditions. The simulation used the 
following cosmological parameters: h — 0.73, Q m = 0.25, 
Q v = 0.75, tt b = 0.045, erg = 0.9, n = 1. The simulation 
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used force softening on a scale of 5 comoving kpc. This 
causes an artificial smoothing of the cores of dark matter 
halos, and as a result, density measurements within ~ 50 
kpc of a halo's core are unreliable. 

The galaxy c atalogues used in th i s pap er consist of 
lightcones within De Lucia fe Blaizotl l|2007t )'s catalogue 
prepared bv lHilbert et al.l (|2009h . These catalogues are com- 
plete for M s teiiar > 10 9 Mq and consist of thirty-two 16 deg. 2 
fields. We split these catalogues into sub-catalogues of "lens" 
and "source" galaxies. The lens catalogues consist of all 
galaxies in the parent catalogue with 0.05 < z < 0.8, and 
the source catalogues consist of all galaxies in the parent 
catalogue with r < 24.5, approx imately the lim i t for shape 
data in the CHFLenS survey mans et alj|2012h . We 
use all galaxies in the lens catalogue for our ray-tracing 
simulations, but only those galaxies with r < 19.8 for 
analysis, m atching the comple teness limit for the GAMA- 
II survey (|Driver et al l l201lh . The source galaxies have 
sh ears calculated from a ray-tracing simulation performed 
bv lHilbert et all which differs from ours in that it uses the 
positions of dark matter particles within the Millennium 
Simulation rather than our halo models (described in Sec- 
tion [2]2]). As such, it provides a useful comparison to check 
that our models generate reasonable results. 

The resolution of redshift and mass within these cata- 
logues was limited by the number of snapshots of the origi- 
nal simulation which were saved. Sixty-four snapshots were 
saved in total, of which 23 are below redshift 1. 

The dark matter halos of the Millennium Simulation 
were populat ed with galaxies through a semi-analytic model 
developed bv lDe Lucia fc Blaizotl |2007l 'l. Dark matter halos 
are identified through a friends-of-friends algorithm, and suf- 
ficiently massive halos are populated with gal axies. Halos are 
also a nalyzed with the SUBFIND algorithm jSpringel et al.l 
l200ll ) to identify subhalos, and self-bound subhalos are sim- 
ilarly populated with galaxies. All of these galaxies are la- 
beled with the ID of the dark matter halo or subhalo that 
contains them. Halos in adjacent snapshots are identified 
as descendants or progenitors if they share at least 50% of 
their dark matter particles, and the galaxies they contain 
are similarly linked. 

We identified galaxies in these catalogues as either "cen- 
tral," "satellite," or "field" galaxies. Central galaxies are 
galaxies which share a h ost halo with at least on e other 
galaxy and are labeled bv lDe Lucia fc Blaizotl (|2007f ) as be- 
ing at the centre of the group's dark matter halo. Satellite 
galaxies are galaxies which share a host halo with at least 
one other galaxy and aren't labeled as central galaxies. Field 
galaxies are galaxies which do not share a host halo with any 
other galaxy: they are centrals in a halo with no satellites. 



2.2 Models 

2.2.1 Halo Models 



distributio ns of all lens halos as truncated NFW profiles, as 
defined bv lBaltz et alj ((20091, with a density profile: 



p(x) = 



Mo 



47rr| x(l + x) 



(1) 



where M = A/200 X (ln(l + c) - c/(l + c)) \ r„ = r 200 /c, 
and x = r/r s . c is the concentration parameter, in princi- 
ple unique to each halo, and r = r t id a i/^s determines the 
truncation radius. This model provides an analytic form fo r 
the weak lensing signal, as outlined in iBaltz et alj (2009). 
In this paper, we assume that for halos that haven't been 
tidally disrupted, r = 2c, which implies r t dai = 2r2oo . This 
assum ptio n is consistent with resul ts from Hilbert fc White] 
lj2010h and lOeuri fc Hamanal (|201lh . the latter of which sup- 
ports a value of r between 2c and 3c. We deter mine c using 
the m ass-concentration relation determined by iNeto et al.l 
l|2007l ) for halos in the Millennium Simulation: 



4.67 x 



M 2 



VQ 1A h-^Uc. 



(2) 



To see how weak lensing signals depend on the mass 
distribution within groups, we experimented with different 
methods of apportioning a group's mass to its constituent 
galaxies. |Pe Lucia fc Blaizotl |2007l Vs catalogue provides the 
total mass of a group's halo (which includes the mass of all 
subhalos) at any snapshot and the stellar masses of all galax- 
ies. Note that the stellar mass is not expected to evolve sis 



nificantly after a galaxy joins a group |Padilla et alj 2009 
so if there is a monotonic relation between infalling halo 
mass and stellar mass, the latter can be used to estimate 
the mass of the galaxy's dark matter halo prior to join- 
ing the group. Here we make this identification using the 
method of "abundance matching" , o r iginally applied t o en- 
tire halos llMarinoni fc Hudsonll2002l: l Yang et al.ll2003l) an d 
later subhalos (|Vale fc Ostrikej 120041: IConrov et al.l |2006| >. 
We use the following formula, from lGuo et al.l (|2010h : 



0.129 x 



fTlstcllar 



rahalo 

1O U - 4 M 



Wi-halo 
10 114 M© 



(3) 



See Appendix I All for details on how we used this formula 
to calculate Mhaio given A/ stc ]iar0 From here on, we will be 
referring to the halo mass calculated this way as the "infall 
mass." The mass of any given group is then the total mass 
of all of its constituent galaxies' halos: 



M to 



Em ha 



(4) 



This prescription, which is based on estimates of halo 
mass from stellar mass, tends to overestimate the to- 
tal group mass, as compared with the group masses 



For ease of calculations, we assume that all dark matter lies 
within spherically-symmetric halof]. We modeled the mass 



Triaxial models ll Jine fc Sutol l2002h provide a more accurate 
representation of real dark matter halos, but this is not necessary 
for our purposes. Per equation lilt , with a sufficient number of 
lenses and sources stacked together, the lensing signal of a triaxial 



halo (which appears elliptical in projection) will be indistinguish- 
able from the lensing signal of a spherically-symmetric halo. This 
is due to the fact that the noise due to neglecting triaxiality is of 
order ~ 20% per pair, while shape noise is 1-2 orders of magni- 
tude larger. 

2 There is a log- normal scatter in this relation of ~ 0.17 dex 
llYane et alj l2012f) , which we do not simulate. Testing showed 
that the effect of this scatter on our results was negligible. 
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in De Lucia fc Blaizotl (|2007T )'s catalogue, as iGuo et al.l 
(|2010l Vs formula was determined based on observational 
data rather than data from the Millennium Simulation. 



2.2.2 Stripping Models 

Using the data for total group mass and satellite infall mass, 
we constructed the following three models for the mass dis- 
tributions of satellite and group halos: 

• Pure Group: A test case in which all group mass is 
assigned to the central halo, and the masses of all satellite 
halos are set to zero. 

• No Stripping: A model expected to correspond to unre- 
laxed groups. Satellites retain their infall masses. The cen- 
tral halo's mass is then set equal to M to t minus the mass in 
satellite halos. 

• Stripping: A model expected to correspond to relaxed 
groups, on average. Satellite halos' masses are decreased rel- 
ative to their infall mass by an amount dependent on their 
projected distance from the group centre. The central halo's 
mass is then set equal to Mtot minus the mass that remains 
in satellite halos. 

We also use the following two datasets for comparison 
purposes: 

• Matched Field: A selection of field galaxies, matched in 
redshift and stellar mass to the satellite galaxies used in the 
above three models. This dataset allows us to measure what 
the lensing signal around the galaxies in the No Stripping 
model would look like if the contribution from being in a 
group environment were removed. 

• Particle Ray- Tracing (PRT ): This dataset uses the ray- 
tracing simulation performed bv lHilbert et~aH (|2009l V which 
uses the positions of dark matter particles rather than our 
ha lo models. Groups are sel ected using the group masses 
in |Pe Lucia fc Blaizotl |2007u 's catalogue, rather than our 
calculated M to t- 

In order to mode l stripping within groups, we applied 
results from Fig 15 of iGao et al. (|2004l ) to estimate a satel- 
lite's retained mass fraction from the projected distance be- 
tween it and its group's centre, using the following function, 
fit by hand to the figure: 



-0.464 



(— y 

V r 2 oo / 



1.03 f-L- 

Vr 2 oo 



+ 0.058. 



(5) 



Minit V r 2 oo - 

We modeled the reduction of satellite mass due to strip- 
ping by decreasing their tidal radii. The relationship between 
mass within the virial radius and tidal radius is given by: 



M ha l = M - 



[(r 2 - l) ln(r) +T7T- (r 2 + l)] . (6) 



As this equation is not analytically invertible, we use a 
solution-space search to find the value of r which gives the 
proper value for 



Minit 



As has been shown by IGao et alj (|2004h . the fraction 
of group mass contained within satellites at a given radius 
increases with distance from the group centre. This implies 
that in the Stripping and No Stripping models, we should 
use a halo density profile for the central halo which con- 
verges to the original profile near the core, but is progres- 
sively lower than the original profile as the radius increases. 




r (kpc) 



Figure 1. Top panel: An illustration of our "contracted NFW" 
profile (red dashed) compared to the initial profile (black solid), 
with a group of Mtot = 10 13 M Q and M sat /M t ot = 0.25. The 
profile converges with the initial NFW profile near the core, but 
falls below it away from the core. Bottom panel: The resultant 
fraction of mass assumed to be in satellites at a distance r from 
the center of this profile. This curve meets our requirements of 
decreasing to zero at the core and rising with radius. 



We accomplish this with the following, modified version of 
the truncated NFW profile, which we call the "contracted" 
NFW profile: 



p{x) = 



M' fr 



47rr^ 3 x'{l + x') 



(7) 



where f r is the fraction of mass retained by the group halo, 
and Mq, r' s , x', and r' are calculated with new concentration 
c' > c, which satisfies the equation: 



frC K 



ln(l + c) - c/(l + c) ln(l + d) - c'/(l + d) ' 



(8) 



This ensures that in the region near the core, where 
i < 1, this profile will converge with the profile given by 
equation (0. Since this profile has a higher concentration 
than the original profile, the density will decrease faster with 
radius than in the original profile. There is no strong physical 
evidence for this specific model; it is used simply because 
it meets the basic requirements of converging to the initial 
NFW profile near the core while having a greater fraction of 
mass in satellites far from the core, and because its lensing 
signal can be analytically calculated in the same manner as 
an NFW profile's lensing signal. 

An illustration of this transformation, along with a plot 
of the mass fraction assumed to be contained within satel- 
lites under this model, can be seen in Fig. [T] The weak lens- 
ing signal for this model can be calculated similarly to the 
original model, simply by using a different concentration and 
applying f r as a scaling factor. For consistency, we use this 
profile for all halos, with f r — 1 for all satellite and field 
halos. 
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3 LENSING SIMULATIONS 

In this section we discuss how we simulate shear from weak 
gravitational lensing and our methods to adjust the lensing 
signals around satellite galaxies to remove the contributions 
from group halos. In Section[3]T] we provide a brief overview 
of weak gravitational lensing and some useful formulae used 
for it. In Section T3. 2 1 we explain our algorithm to determine 
the ellipticities of source galaxies. In Section T3. 31 we explain 
how we analyze the data to calculate lensing signals, and we 
discuss both of our methods for correcting satellites' lensing 
signals for the contribution from group halos. 



3.1 Lensing Overview 

Here we briefly review weak lensing. A more thorough ex- 
planation of the tools of weak lensing , with an emphasis 
on len sing in clusters, can be found in iKneib fc Nataraianl 
(2011), but a brief explanation of the physics and terminol- 
ogy follows: the presence of mass distorts spacetime, causing 
light to be deflected when it passes near a massive object. 
In the so-called "weak lensing" regime, this manifests as a 
coherent distortion of the shapes of background galaxies, 
stretching them tangential to the foreground object, in an 
effect known as "shear." Although the variation in the in- 
trinsic shapes of background galaxies is significantly larger 
than the shear caused by foreground objects, the shape of 
any given background galaxy will still provide an unbiased, 
albeit low signal-to-noise, estimate of the local shear. This 
makes it possible to gather usable data by stacking a suffi- 
cient number of lens galaxies together. 

In practice, shear is usually observed and measured 
through the dimensionless parameter 7t (known as the tan- 
gential shear), but as the magnitude of this parameter de- 
pends on the redshifts of both the source galaxy and the 
foreground object, for this paper we will use the measure- 
ment AE, which is related to jt through: 



<7t> 



AS 

Scrit 



in a circle of any size around any given point, where: 
c 2 A 



(9) 



(10) 



4ttGAsA ' 

and D B is the angular diameter distance to the source, A 
is the angular diameter distance to the lens, and As is the 
angular diameter distance from the lens to the source. 

For an arbitrary mass distribution, (AE) at a radius r 
from a central point can be calculated through: 



(AE(r)) = E(< r) — E(r), 



(11) 



where E(< r) is the surface density averaged for all points 
contained within this circle, and E(r) is the surface den- 
sity averaged for all points on the edge of this circle. This 
prescription works even for mass distributions that are not 
axisymmetric, as long as all points in a given annuli around 
a lens object are stacked together. 



3.2 Determining Ellipticity 

While full ray-tr acing, such as that performed by 
iHilbert et af] (2009), is a powerful tool for estimating weak 



lensing signals caused by an arbitrary mass distribution, it 
has a few limitations. The accuracy of the ray-tracing is 
limited by the resolution of the simulations used to generate 
the lens mass distribution, which causes the algorithm to 
noticeably underestimate the lensing signal within 50 kpc 
of halo cores (this effect can be seen in Fig. [3}. Addition- 
ally, the algorithm has no method to easily model different 
hypothetical mass distributions, as we wish to test here. As 
such, we've developed a modified ray-tracing algorithm, in 
which we assume that all dark matter lies in spherical halos 
in amounts and distributions dependant on the models de- 
scribed in Section [2.2.21 The algorithm proceeds as follows: 

(i) Lens galaxies are assigned mass, concentration, tidal 
radius, and scaling factor under one of the models described 
in Section |2~2~21 

(ii) In order to eliminate edge effects, the catalogue of lens 
galaxies used for this algorithm includes galaxies within a 
wider field of view than the catalogue of source galaxies. 
Since the fields are 4x4 degrees, we accomplish this by 
only using source galaxies within a central 1.4 x 1.4 degree 
square field, which also significantly reduces the amount of 
computational time needed for later analysis. For the rest of 
our analysis, we also only use lens galaxies that lie within 
this central field. 

(iii) Each source galaxy is initialized with zero ellipticity 
in both components (shape noise is simulated at a later stage 
of analysis). 

(iv) For each source galaxy, tangential 
shear is applied to it for every lens galaxy 
where ([RAi cns — RA sourcc ] x cos(Deci ons )) 2 + 
(Deciens - Dec sourcc ) < (A max ) and 

^lcns ^ 2-sourcc - 

A max is an upper limit to conserve computational time and 
ensure that the rectangular boundaries of the field will not 
cause artifacts in the lensing sig nal. The str e ngth of the 
shear applied is determined using Baltz et al. (2009)'s cal- 
culations for a truncated NFW halo and our modifications 
for a contracted NFW profile (see Section [2.2.21) . 

(v) Shear is added linearly to the source galaxy's elliptic- 
ity for all lens galaxies within this lightcone. 

This process in principle accounts for the two-halo term 
seen in the lensing signal around groups, which is caused by 
other nearby groups. However, our application of a cut-off 
angle A max suppresses lensing signals at large radii, mean- 
ing the two-halo term is not observable in our simulated 
lensing signals. This method does not account for dark mat- 
ter that is correlated with the positions of galaxies and 
groups but not considered part of their halos, which also 
results in a slight suppression of the lensing signal at large 
(R J> 1000 kpc) scales. Since we are interested in the strip- 
ping of subhalos, for which the relative signal is strongest at 
intermediate (50 kpc < R < 400 kpc) scales, this suppres- 
sion at large scales will have no effect on our analysis. 

3.3 Stacking Lensing Signals 

Once the source catalogue has been prepared with simulated 
ellipticities for each source, we proceed to stack the lensing 
signal around a sample of lenses in order to observe the 
average signal for the sample. We separate our sample of 
lens galaxies into samples of central galaxies and samples of 
satellite galaxies. 
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3.3.2 Ensemble Method 

The second method involves estimating the masses of groups 
and their contributions to the lensing signals around their 
satellites. First, we create a sample of groups with M m i n < 
M < Af max and each with at least one satellite with m m - m < 
m < m max . We use a weighted average of the groups' masses 
to generate an "ensemble" halo, representative of the sam- 
ple. We then assign this ensemble halo a concentration based 
on equation (5). For each source in the sample, we then 
use our ensemble halos and the positions group centres to 
calculate the group's expected contribution to the source's 
ellipticity and subtract off this amount. 



Figure 2. An illustration of the Mirror Method. For each satel- 
lite used as a lens, we also calculate the lensing signal around a 
point opposite the group centre from this satellite. We subtract 
this signal from the signal around the satellite to correct for the 
contributio n of the group's halo to the lensing signal. See also 
Figure 1 o f IPastor Mira et al.l |201lh and related discussion. 

With satellite galaxies, we are interested in comparing 
the lensing signals of similar satellites in groups of different 
masses. For the sake of clarity in discussing these limits, 
we use m m in and m m ax to define the lower and upper limits, 
respectively, for satellite mass in our samples, and M m i n and 
M max to define the lower and upper limits, respectively, for 
the masses of the groups in which they reside. 

With each subset of lens galaxies, we calculate the aver- 
age tangential shear within radial bins. The results of these 
calculations and analysis of the results can be seen in Sec- 
tion [472] and Section [4. 31 The subsets of satellites will typi- 
cally show some contribution from their groups' halos to the 
lensing signals around them. Since the observed lensing sig- 
nal around satellites is the combination of the contributions 
from the satellites' halos and their groups' halos: 

AEobs = ASgroup + AE sa t, (12) 

we can isolate the portion of the signal contributed by the 
satellites' halos (AE sa t) by estimating the groups' contri- 
bution to the lensing signal (AE group ) and subtracting it 
from the observed lensing signal (AE b s ). We have applied 
the two methods detailed below to do this, and we compare 
their results in Section |4.3. f I and Section |4. 3. 21 

3.3.1 Mirror Method 

The first method for isolating the satellite's contribution 
to the lensing signal involves using a "mirror" point in the 
group to measure the group's contribution to the satellite's 
lensing signal. This sample point is placed opposite the 
group centre from the satellite, at the same distance from it, 
as illustrated in Fig. [2] In an ideal scenario, the lensing sig- 
nal around this sample point will include exactly the same 
contribution from the group's halo as at the satellite's lo- 
cation, while including little contribution from the satellite 
itself. 

This method is in practice identical to the "Galaxy- 
Galaxy Lensing with Calibration" method used by 
IPastor Mira et all (|201 lh . Figures 1 and 2 in that paper pro- 
vide a useful illustration of how this method functions, and 
discussion of it is found in Section 3. 



4 LENSING ANALYSIS 

In this section we present the lensing signals that result from 
our simulations and discuss the effects of various sources of 
error on these signals. In Section [4.11 we discuss the domi- 
nant sources of error in this type of analysis. In Section |4. 21 
we present the lensing signals that result when group centres 
are stacked together. In Section [4.31 we present the lensing 
signals that result when satellite galaxies are stacked to- 
gether, we show the effects of our methods to reduce the 
contributions of group centres to these signals, and we in- 
vestigate the effects of centering and group mass errors on 
these signals. 



4.1 Dominant Sources of Error 

In order to properly assess whether our methods for measur- 
ing satellites' lensing signals are viable with observational 
data, it is necessary to account for the sources of error 
present in such datasets. In this section, we explain the 
major sources of error we expect in datasets and how we 
replicate them in our simulations. 



4-1.1 Shape Noise 

There is significant variation in the orientation of galax- 
ies in the sky, independent of the effects of gravitational 
lensing. The distribution of the unlensed shapes of back- 
ground galaxies can be approximately modeled by a Gaus- 
sian distribution for both components of ellipticity, with a 
mean of zero and a standard deviation of 0.28. This value 
was calculated from the distribution of ellipticities measured 
in background gala xies in the CFHTLenS, estimated with 
the lensfit method dMiller et all 120071 ; iKitching etHI 120081 ; 
iHevmans et alj|2012h .~ 

We account for the impact of shape noise on our mea- 
surements by setting the size of the error bars in our plots 
to account for its expected impact. This allows us to pre- 
cisely predict the mean expected measurements, unaffected 
by shape noise, and to estimate the scatter that will be seen 
in observational data. Our error estimates assume data from 
a survey similar to the overlap of the CFHTLenS survey 
and the GAMA-II survey. This hypothetical survey covers 
50 deg. 2 of the sky, with an effective 15 sources per arcmin. 2 
and galaxy spectra complete to r < 19.8. 
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4-1-2 Group Mass and Group Assignment Errors 

With observational data, group masses must be estimated 
through one of various methods, such as using the veloc- 
ity dispersion of constituent galaxies, the luminosity of the 
group, X-ray emissions, or lensing signals. In the GAMA-II 
survey, group masses will be estimated primarily through 
velocity dispersions, so this is the method that we simulate 
here. To do this, we first determine which of a group's galax- 
ies are within the survey's detection limits. For each of these 
galaxies, we generate a random line-of-sight velocity from a 
Gaussian distribution, where the disp ersion is determined 
by th e group's total mass, as given in iDe Lucia fc Blaizotl 
l|2007|y s catalogue. For all calculations relating ai and M200, 
we assume a Singular Isothermal Sphere (SIS) profile, which 
gives us the relation: 

o-i = (lOG/iMaoo) 173 /^- (13) 

When groups are detected through a friends-of-friends 
algorithm, it is expected that some portion of the galaxies 
believed to be within a group will in fact be "interlopers." 
This happens when field galaxies lie at a similar position in 
the sky as a group at a different cosmological redshift, but 
with a peculiar velocity that makes their apparent redshift 
lie within the redshift distribution of galaxies within the 
group. For a given FoF-detected group, typically 20% of the 
galaxies identified as b elonging to it are in fact interlopers 
[ Robotham et alJl201lh . 

To simulate the effects of interlopers on group mass es- 
timates, we first determine the number of interlopers that a 
given group contains. Since the number of interlopers should 
be proportional to the projected area of the group, we set the 
mean number of interlopers in a given group proportional 
to A^ 3 . We then use a Poisson distribution to generate the 
number of interlopers for each group, normalized so that a 
group with 10 members will have an expected 2.5 interlopers 
(resulting in a 20% interloper fraction). For each of these in- 
terlopers, we then generate a velocity for them, drawn from 
a uniform distribution, u min — v t ol < "rand < « mK + Wtoi, 
where u m i n and v max are the minimum and maximum veloc- 
ities of real members of this group, and v to i is the velocity 
tolerance used by the FoF algorithm. In the algorithm used 
for the GAMA-II survey, v%oi is riot fixed, and in fact de- 
pends on the luminosities of the galaxies being linked, but 
for simplicity's sake, we use a constant value of 150 km/s 
here. 

Once velocities for real members and interlopers are 
generated, we then calculate the predicted group mass using 
these velocities. When this mass is used in our analysis, we 
use the subscript "dyn." 

The presence of interlopers will also affect the lensing 
signal for the "Stripping" model, as interlopers would typi- 
cally be field galaxies, which will exhibit an unstripped mass 
distribution. To simulate this effect, we recalculate the signal 
in the "Stripping" model as: 

AEs.new = 0.8 x AEs.om + 0.2 x AE NS , (14) 

where AEs is the lensing signal found with the "Stripping" 
model, and AEns is the lensing signal found with the "No 
Stripping" model. 

We use both the modifications to group mass, which 
affects the binning of centres, and the simulation of the con- 



tamination caused by interlopers only for the plots shown 
in Section 14.3.21 which attempt to predict the strength of 
the signal with a CFHTLenS + GAMA-II-like survey (see 
Section liTTTT) . 

4-1.3 Centering Errors 

It is not trivial to determine the position of a group's cen- 
tre. Since the majority of the mass in a group is dark mat- 
ter, and galaxies are an imperfect tracer of the dark matter 
distribution, any method which uses galaxies' positions to 
estimate the centre of a group's dark matter halo will have 
some degree of error. The amount of error can be roughly 
estimated by analysis of t he velocity di s persio n about sus- 
pected centres, as done bv lSkibba etafl l|201ll ); by compar- 
ing the results of multiple independant methods for identify- 
ing the group's centre, or by fitting the lensing signal around 
a s tack of groups with a profile convolved with error, as done 
by iGeorge et al.l l|2012T ); or through direct analys is of N- 
body simulations, as done by I Johnston et akl (120071 ). While 
Ijohnston et al.l (|2007l )'s analysis of simulations implies that 
rich groups are less likely to be miscentered than poor 
groups, this d oesn't seem to be the case with real groups. Re - 
cent work by ISkibba et ail (|201ll ) and IGeorge et ail (I2012T ) 
all show that this trend is in fact reversed in real groups. For 
the most massive groups, the centre is misidentified roughly 
40% of the time, impl ying that mis i dentif ying the BCG is 
a significant problem. [George et all (|2012i )'s work has also 
demonstrated that using either the most massive or bright- 
est galaxy in a group as the centre is the optimal method 
when X-ray observations aren't available. 

Therefore, to simulate the effects of centering errors, we 
assume that all groups are miscentered according to a 2D 
Gaussian with a proportional to their scale radii r s . When 
groups are stacked together, this leads to roughly the same 
net effect as observed previously, but it allows us to adjust 
the magnitude of the effect by altering only the constant of 
proportionality (rather than both the fraction of groups in 
which the centre is misidentified and the dispersion of the 
false centres relative to the real centres), to easily compare 
different amounts of error. Th is does not, however, perfectly 
match the best fit models of IGeorge et al.l (|2012h , so only 
qualitative conclusions should be made from our simulations 
involving centering error. We compare models with centering 
errors a — 0, r s , and 2r s , which covers the expected range of 
error. When real data are available, it will be necessary to 
determine the proper amount of error through an analysis of 
the lensing signals around suspected group centres, and then 
recalibrate our pred ictions for it, using a method similar to 
IGeorge et all l|2012i Vs method of fitting the lensing signals 
with an NFW profile convolved with centering error. 

For the Mirror Method, centering errors will have an 
impact because if there's error in the position of the group's 
centre, there will also be error in the position of the mirror 
point. On average, this error will cause the mirror point 
to lie farther away from the group centre than it actually 
should, and so it underestimates the group's contribution to 
the lensing signal around its satellite. 

For the Ensemble Method, centering errors will affect 
the calculated signal from the ensemble halo. The reason for 
this is that errors in the positions of group centres will sup- 
press the group signal at low radii (R <. <7 con ), thus causing 
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Figure 3. Lensing signal, measured with AS, around a sample of 
galaxy groups with 10 13 Mq < Mtot < 10 14 Mq drawn from the 
Millennium Simulation, with no centering or group mass errors 
applied. All four mass models detailed in Section 12.21 are shown 
here, along with the lensing signal that results from using the 
particle ray-tracing (PRT) shear data. The solid line shows the 
predicted shear for a group which has the weighted average mass 
of the sample (Aftot = 4.00 X 10 13 Mq) at the weighted average 
redshift of the sample (z gr = 0.18). Error bars are shown for 
the No-Stripping model which represent the projected errors for 
data from a hypothetical CFHTLenS + GAMA-II-like survey (see 
Section HXTJ. 

a lower mass to be fit. They will also change the calculated 
satellite distances from group centres, increasing them more 
often than decreasing them, and thus decreasing the ensem- 
ble halo's typical impact on source ellipticity. In practice, 
this effect is somewhat smaller than the effect of centering 
errors on the Mirror Method (see Section T4. 3. 21 for details). 



Figure 4. Lensing signal, measured with AS, around a sample of 
satellite galaxies with 10 10 Mq < m^aio (as determined with the 
"No Stripping" mass model), and 10 13 M < M to t < 10 14 Mq, 
drawn from the Millennium Simulation. All three mass models 
detailed in Section l2.2l are shown here, along with the "Matched 
Field" dataset and the lensing signal that results from using 
the particle ray-tracing shear data with a similarly-cut sam- 
ple of lenses. Error bars are shown for the No-Stripping model 
which represent the projected errors for data from a hypothetical 
CFHTLenS + GAMA-II-like survey (see Section HXTJ . 

implies that our contracted NFW profile is not a perfect 
model for the background halo of a group with significant 
substructure in it. More work will be required to determine 
a better model for this. This effect is minor, however, and 
all lensing signals are consistent with an NFW profile. If 
centering errors (see Section I4.1.3P are present in observa- 
tional data, the lensing signal around group centres should 
be suppressed at low radii. 



4.2 Lensing Signal around Group Centres 

Since the expected form of the lensing signal around group- 
centres is well-understood, measuring this signal in our sim- 
ulations provides a useful test to confirm that our models 
behave as expected. Fig. shows the average lensing signals 
around a sample of group centres with 10 13 M < Mtot < 
10 14 Mq, where Mtot is the total mass of the group, for 
all of our mass models and the particle ray-tracing data, 
with no centering or group mass errors applied. Our models 
show rough agreement with the PRT lensing signal at large 
radii, although the PRT signal drops significantly below the 
models' si gnals for sources wit hin » 50 kpc. This effect was 
noticed by Hil bert et al.l (|2009l ) , and is due to force softening 
used within the Millennium Simulation, which flattens the 
cusps of halos, as well as the smoothing applied during the 
ray-tracing to suppress particle shot noise. 

All of our models appear similar in this regime, with 
a slight tendency for the "No Stripping" and "Stripping" 
models to fall slightly below the "Pure Group" model. This 



4.3 Lensing Signal around Satellites 

Fig. |3] shows the average lensing signals around a sample of 
satellite galaxies with 10 10 Mq < mhaio, where rrihaio is as 
used in equation ©, and 10 13 M < Mtot < 10 14 Mq, for 
all of our mass models and the particle ray-tracing data. 
The differences between our models are in principle de- 
tectable with the amount of data available in our hypothet- 
ical CFHTLenS + GAMA-II-like survey (see Section gXTJ. 
At radii larger than ~ 50 kpc, the PRT signal lies close to 
our No Stripping model, which is an encouraging sign that 
our models are giving reasonable results. However, our "No 
Stripping" model is noticeably different from the "Matched 
Field" dataset, due to the contributions of the satellites' host 
groups to the lensing signal, and this difference is greater 
than the difference between it and the "Stripping" model. 
The "Pure Group" model, as shown here, shows only the 
lensing signal that results from setting the satellite masses 
to zero in our simulations, and so it can be used as an ap- 
proximation of the group halos' contribution to the lensing 
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signal around satellites. In order to use the "Matched Field" 
dataset as a standard for comparison, it will be necessary to 
eliminate the groups' contributions to the lensing signals 
measured here. This may be possible through one of the 
methods outlined in Section 13.31 Below we have presented 
the results of our tests of these methods. 

4-3.1 Idealized Group Signal Subtraction 

As detailed in Section ^. 3. II the Mirror Method uses a sam- 
ple "mirror" point opposite the group centre from a satel- 
lite to estimate the group's contribution to the satellite's 
lensing signal. Fig. [S] shows the result of applying this to 
a sample of satellite galaxies with 10 10 Mq < mhaio and 
10 13 M < Mtot < 10 14 M , in the idealized scenario in 
which there are no errors in the positions of group centres 
or the estimated masses of groups. As can be seen in the 
figure, the curve for the "No Stripping" model is very sim- 
ilar to that of the "Matched Field" dataset, and the "Pure 
Group" model has disappeared off the bottom of the plot 
(it is consistent with zero). This implies that the method 
does a satisfactory job of eliminating the groups' contribu- 
tion to the lensing signal around their satellites. The signal- 
to-noise in the difference between the "No Stripping" and 
"Stripping" models is quite strong. 

The Ensemble method, explained more thoroughly in 
Section 13.3.21 uses the average mass of the groups in the 
sample as an average "Ensemble" halo. For each group cen- 
tre, its predicted shear is then subtracted from the elliptic- 
ities of all nearby source galaxies. After this, satellites are 
stacked together, and the lensing signal around them, using 
the altered ellipticities, is measured. Fig. [5] shows the result 
of applying this method to a sample of satellite galaxies with 
10 10 M < m ha io and 10 13 M < M tot < 10 14 M , also in 
the idealized scenario of no centering or group mass errors. 
The signals from the "No Stripping" model show good agree- 
ment with the "Matched Field" model. Notably, the differ- 
ence between the "Stripping" and "No Stripping" models 
here is greater than when the Mirror Method is used. This 
is because the Mirror Method results in some of the lensing 
signal around satellites being lost, particularly for satellites 
near the centres of groups. Although the Ensemble Method 
requires more assumptions about the physical properties of 
the groups, it does not face this problem, and so the resul- 
tant signal-to-noise in the difference between the two models 
is larger. 

4-3.2 Realistic Group Signal Subtraction 

With observational data, the positions of group centres will 
not be precisely defined, so it is necessary to assess the im- 
pact of centering errors on the utility of this method. Figure 
[6] shows how the lensing signal appears if centering errors of 
cr ccn = 2r s are assumed. We simulated this by applying a 2- 
dimensional Gaussian deviate to the positions of groups cen- 
tres. This amount of cent ering error is signifi cantly higher 
than the errors found by I George et al.l (|2012l ). which were 
typically of order <T CC n ~ 0.5r s — lr s , and so this provides an 
upper bound to the amount of centering error that might be 
seen in observational data. As can be seen through compar- 
ison with Fig. O centering errors have a significant impact 
on the resultant lensing signal. 



To determine whether or not we will be able to detect 
a difference between the "Stripping" and "No Stripping" 
models when centering errors are present, we've plotted the 
best-fit NFW mass (determined by minimizing the x 2 statis- 
tic for the group-subtracted lensing signal around satellites 
in radial bins between 10 and 400 kpc) for our "Stripping" 
and "No Stripping" models in both group mass bins in Fig- 
ure [7] for varying amounts of centering error. As this figure 
shows, centering errors increase the best-fit masses in all 
mass bins, but it is only a significant problem for the more 
massive groups. If we have no knowledge of how large center- 
ing errors in observational data will be, the uncertainty they 
cause will be too large for us to discern between stripped and 
unstripped models. This is unlikely to be the case, however, 
as we will be able to measure the lensing signals around 
group centres to estimate the magnitu de of centering error s 
through a fit such as that performed bv lGeorge et alj (2012), 
and significant amounts of error will manifest as a much- 
poorer x 2 fit- Details for all models, assuming no centering 
error, can be seen in Table [1] 

Figure [S] shows the projected best-fit NFW masses for 
different group mass bins if there is significant error in our es- 
timates of group masses. This causes bleeding between bins, 
increasing the expected mass in all group mass bins. We've 
illustrated this effect by using the dynamical mass for bin- 
ning. Note that estimates of halo mass from stellar mass will 
likely have much lower scatter, so this should be seen as an 
upper bound for the effect of binning error. 



5 DISCUSSION 

As our tests have shown, tidal stripping will have a signifi- 
cant effect on the lensing signals of satellites within groups. 
However, various sources of error make the detection of this 
effect more difficult in practice. As with any weak lensing 
measurements done, shape noise is the largest source of error 
and the most difficult to overcome - it can only be reduced 
by gathering more data. By determining the expected sig- 
nals for stripped and unstripped satellites, we can estimate 
the amount of data that will be necessary to discriminate 
between stripping models. 

To measure tidal stripping, it's necessary to eliminate 
the contribution of the shared group halo to the lensing sig- 
nal around satellites. Any method to do this will introduce 
its own sources of error, which can be comparable in magni- 
tude to shape noise. In the cluster regime, it's possible to get 
very accurate mass estimates for individual clusters, and so 
these errors can be minimized even for individual clusters. In 
the group regime, it is only possible to estimate group mass 
sufficiently well when multiple groups are stacked together, 
but the far greater numbers of groups make this feasible. 

5.1 Optimal Methodology 

Our tests of both the Mirror Method and Ensemble Method 
show that in ideal circumstances, they both successfully 
eliminate the groups' contribution to the lensing signals 
around satellites. In this scenario, the Mirror Method is less 
biased, but suffers from a greater error due to shape noise 
- a result of the fact that it relies in the difference in signal 
between two points, causing the error to be increased by a 
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Figure 5. Same as Fig. [4] except with the Mirror Method (left panel) and Ensemble Method (right panel) applied to reduce contribution 
to the lensing signal from groups' halos (except with the "Matched Field" dataset). 
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Figure 6. Same as Fig. [4] except with the Mirror Method (left panel) and Ensemble Method (right panel) applied to reduce contribution 
to the lensing signal from groups' halos (except with the "Matched Field" dataset), assuming centering errors of a = 2r s . 



Table 1. Statistics for our galaxy samples and NFW profile fits, using the Ensemble Method and assuming o" CC ii = 0. All masses are in 
units of 10 10 Mq. S/Np_s is the signal-to-noise ratio in the difference between the fitted masses for the "Field" and "Stripping" models, 
and S/Nns— s i s * ne same for the difference between the "No Stripping" and "Stripping" models. 
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Figure 7. Best-fit NFW masses for a sample of satellite galaxies with 10 10 Mq < mi la i (as determined with the "No Stripping" mass 
model) drawn from the Millennium Simulation, using only source galaxies within 10 to 400 kpc of the lens galaxies. Satellites from three 
ranges of group mass are included, shown at the average mass for the range, and the Mirror Method (left panel) and Ensemble Method 
(right panel) have been applied with three different levels of centering error assumed. Centering errors used were <r ccn = (left), cr ccn = r B 
(centre), and <T C en = 2r s (right). Error bars shown are the projected errors for data from a hypothetical CFHTLenS + GAMA-II-like 
survey (see Section 14.1,1 [I . The "Matched Field" mass is observed to increase with group mass as more massive groups typically contain 
more massive satellites, and these more massive satellites are then matched to more massive field galaxies. 
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Figure 8. Same as Fig. [7] except satellites have been binned by their groups' estimated dynamical masses instead of the actual masses. 



factor of v2- The Ensemble Method causes errors to be con- 
tributed from the estimate of the average group mass, but 
the total errors still tend to be less than the errors from the 
Mirror Method. 

When comparing the "No Stripping" and "Stripping" 
models, we see that the difference between the two models is 
larger when the Ensemble Method is used. This is due to the 
fact that the Mirror Method, by subtracting out the signal 



from around sample points, also subtracts out a portion of 
the lensing signal around satellites. In particular, for a point- 
mass satellite a distance R from the centre of its group, the 
Mirror Method will detect zero signal for it on average for 
sources at annuli of radius greater than 2R (the distance 
from the satellite to the sample point). The satellites which 
are found closest to group centres tend to be the ones which 
have undergone the most stripping, but these are also the 
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satellites which lose the largest amount of their own lensing 
signals when the Mirror Method is used. This suppresses 
the difference between the "No Stripping" and "Stripping" 
models. 

Although neither method works perfectly when center- 
ing errors are present, the Ensemble Method shows some- 
what less sensitivity to centering errors than does the Mir- 
ror Method. This difference is most significant for centering 
errors of a — r s . Note that we expect to be abl e to quan- 
tify centering errors through [George et al l l|2012ft 's method 
of fitting the lensing signal around group centres with an 
NFW profile convolved with centering error, but it is still 
useful to minimize the expected offset due to centering er- 
rors. 

Overall, the Ensemble Method is preferable to the Mir- 
ror Method. It presents a better signal-to-noise in the dif- 
ference between the expected signals for stripped and un- 
stripped satellites, and it is less vulnerable to centering er- 
rors than the Mirror Method. 

Other methods, such as a maximum-likelihood analysis, 
which work similarly to the Ensemble Method, will likely 
have a similar signal-to-noise to what we've calculated here. 



5.2 Prospects for Stripping Detections 

A detection of tidal stripping with either of the methods 
we have proposed will first require the issue of centering 
errors to be addressed. If they can be addressed a nd co r- 
rected for (eg. by using the method of lGeorge etafl |2012l )). 
the Ensemble Method then provides the better S/N. By us- 
ing the Ensemble Method and comparing the S/N in the 
difference between the Stripping and No Stripping models, 
we can estimate how much data will be needed in differ- 
ent group mass regimes in order to differentiate between the 
models. Table [1] lists our predicted S/N values for a mea- 
surement using the overlapping region of the CFHTLenS 
and GAMA-II surveys. We list S/N values for both the dif- 
ference between the "No Stripping" and "Stripping" mod- 
els, and the difference between the "Matched Field" dataset 
and the "Stripping" model. Of these two, the latter compar- 
ison is what we will measure in observational data, but it 
may be influenced by systematic errors in the application of 
the Ensemble Method. The former comparison is less likely 
to be influenced by systematic errors, as we used the same 
fields and realization of noise for both the "No Stripping" 
and "Stripping" models. As such, it provides a more realistic 
prediction for what we might see with observational data. 

Of the three mass regimes tested, we see the highest 
S/N, 1.88, in the group mass bin of 10 13 M Q < M to t < 
10 14 Mq . Although this doesn't quite reach the 2a threshold 
for a detection, it is still worth making this measurement for 
various reasons: 

• We do not apply any weighting scheme to our lenses or 
sources. An optimal weighting scheme could serve to boost 
the S/N. 

• When all three mass bins are combined together, the 
predicted S/N rises to 2.86. This will allow us to detect the 
presence or absence of stripping, but not the group mass 
regimes in which it is taking place. 

• If stripping does occur and is stronger than predicted 



by our models (see equation (J5J), a detection will be more 
likely. 

• Even if no 2a detection can be made from the 
CFHTLens+GAMA-II dataset, the results of this measure- 
ment can be combined with future measurements to poten- 
tially add significance to those results. 

Our predicted significance can be improved 
upon if data from other surveys is combined with 
the CFHTLenS+G AMA-II dataset. Of note, KIDS 
jde Jong et al l I2012T I will overlap the GAMA-I survey in 
150 deg. 2 when it is complete, and we estimate it will 
contain an effective 8 sources per arcmin. 2 . Assuming 
that the overlap between KIDS and GAMA-I will be 
80% unmasked, this results in a dataset that is ~ 1.3x 
larger than CFHTLenS+GAMA-II. Table [2] shows the 
predicted significances for the three mass regimes we've 
tested, assuming the KIDS+GAMA-I dataset is used alone, 
and also combining it with the CFHTLenS+GAMA-II 
dataset. KIDS+GAMA-I alone is somewhat better than 
CFHTLenS+GAMA-II, and the combination of the two 
datasets has a > 50% chance of a 2a detection in each mass 
bin, if stripping is occurring at the strength we've modeled. 

In conclusion, we will soon be able to use weak gravi- 
tational lensing to measure tidal stripping within even low- 
mass galaxy groups. A better understanding of tidal strip- 
ping will hopefully allow us insight into its role in the evo- 
lution of galaxies within the group environment. 



6 CONCLUSIONS 

In this paper, we have shown that for analysis of the lens- 
ing signals of satellite galaxies, the Ensemble Method (see 
Section I3.3.2[) successfully subtracts out the contributions 
of these satellites' host groups to their lensing signal. This 
method performs better overall than the Mirror Method (see 
Section I3.3.1|) . 

We have shown here than errors in identifying the cen- 
tres of group halos have a significant impact on the analysis 
of satellites' lensing signals. The contribution of these center- 
ing errors to the lensing signals around satellites is typically 
larger than the effect of tidal stripping on these signals. It is 
thus necessary to properly estimate the amount of centering 
error present in a sample before an assessment can be made 
of the presence or absence of tidal stripping. 

If centering errors in the dataset can be measured, up- 
coming data from the GAMA-II survey, when combined with 
the CFHTLenS dataset, may be able to detect tidal strip- 
ping within groups of mass 10 13 M < M tot < 10 14 M to 
2a significance, depending on the strength of stripping in 
reality and the scatter of the observational data. Data from 
KIDS and GAMA-I can be used to improve the significance 
of this measurement, and to also potentially detect stripping 
in lower-mass groups. 
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Predicted S/N NS _s 



logM 
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CFHTLenS+KIDS+GAMA 
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1.43 


1.63 
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1.88 


2.14 


2.85 
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1.61 


1.84 


2.44 
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APPENDIX A: APPROXIMATIONS USED 
Al Deriving Halo Mass from Stellar Mass 



In Section 12.21 we found it necessary to estimate the mass 
of a galaxy's dark matter halo, given its stellar mass. While 
equation ([3]) can be easily used to estimate a galaxy's stellar 
mass if its halo mass is known, the equation is not invert- 
ible. We accomplish our task of estimating Mh a i through 
an iterative process: 

(i) Begin with an estimate of the galaxy's halo mass, 
Mjnit. For this, we use the M 200 value for the galaxy given 
in |Pe Lucia fc Blaizotl ([2007) 's database. 

(ii) Set M_i, M-2, M_ 3 , and M_ 4 all equal to M ini t. 
(hi) Set Mt cs t equal to (M_i + M_ 2 + M_ 3 + M_ 4 ) /4. 
(iv) Determine M new through: 



M new = 



A^stcllar 

0.129 



+ 



10 11 - 4 M 



10 11 - 4 M 



(Al) 



(v) If M now , M_i, M_ 2 , M_ 3 , and M_ 4 all differ by no 
more than 0.1%, return M ncw as Mhaio. 

(vi) Set M_4 = Ms, Ms = M-2, M-2 = M-i, and 
M_i = M ncw 

(vii) Repeat steps 3 to 6 until M has converged, for a 
maximum of 1,000 loops 

(viii) If M has not converged, return Mi n i t as Mhaio 

The use of multiple past estimates of Mhaio in our code 
was found to be necessary in order to ensure convergence, 
as when fewer past estimates were averaged and used, the 
algorithm would often end up alternating between two vastly 
different estimates for Mhaio- With the algorithm as written, 
Mhaio was found to converge in every case tested in a sample 
4x4 deg. 2 field. 

A2 Altering Concentration of NFW Profiles to 
Preserve Cusp Distribution 

In Section 12.21 we mentioned an alteration we made to 
the NFW mass profile which leaves the density profile un- 
changed near the centre, but is lower than the initial pro- 
file at large radii. To illustrate how this works, examine the 
formula for the NFW density profile, with S c , r 3 , and x ex- 
panded: 



p{x) = 



M 2 oo*(ln(l + c)-e/(l + e))~ 1 



V I"200 / V V r 200 > I 



In the limit r <C r2oo / c, this becomes: 



p{x) = 



M200 1 



47rr 200 r ln (! + C ) _ C/C 1 + C ) ' 



(A2) 



(A3) 
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Figure Al. Percentage error in our estimates of c', as calculated 
from equation IA4t . as a function of lc (defined in equation |(A5j) • 
The error is less than 2% for all reasonable values for the group's 
initial concentration c and the fraction of mass f r retained by 
the group's halo. Significant error is only seen for extremely low 
concentration halos, which aren't modeled in our simulations. 

Since we modify the NFW profile by introducing a scale 
factor f r to the profile and modifying c, equation ((8| must 
be satisfied for the density in the core region to remain un- 
changed. As this equation is not solvable for c', we use the 
following approximation: 



c = 0.1171c 4 - 1.70k 3 + 11.10/c 2 - 30.7k + 32.8, (A4) 



This approximation has less than 2% error for 3 < lc < 
7, as can be seen in Fig. IA21 which covers the range of lc 
possible for 4 < c < 15 and 0.1 < f r < 1. 



where: 




(A5) 



